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Abstract 

The presence of the aging phenomenon in the homogeneous cooling state (HCS) of a granular 
fluid composed of inelastic hard spheres or disks is investigated. As a consequence of the scaling 
property of the A-particle distribution function, it is obtained that the decay of the normalized two- 
time correlation functions slows down as the time elapsed since the beginning of the measurement 
increases. This result is confirmed by molecular dynamics simulations for the particular case of 
the total energy of the system. The agreement is also quantitative in the low density limit, for 
which an explicit analytical form of the time correlation function has been derived. Moreover, 
the reported results provide support for the existence of the HCS as a solution of the N-particle 
Liouville equation. 

PACS numbers: 45.70.-n,51.10.+y,05.20.Dd 
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I. INTRODUCTION 



In many far from equilibrium states, it has been observed that the relaxation (or response) 
rate decreases as the "age" of the process increases. Here "age" refers to the time elapsed 
since the beginning of the considered experiment. Then, it is said that the system ages, or 
that it exhibits an aging phenomenon. A revision of the concept of aging in spin glasses and 
in other systems can be found in ref. 

One of the typical experiments for the analysis of aging in a given physical system is the 
study of the two-time correlation functions of some of its properties, since they are related 
with the response of the system to a given perturbation. Let Cab^w, t) denote the two-time 
correlation function of the magnitudes A and B of the system, the former being measured at 
time t w and the latter at time t >t w . In a system at equilibrium, CAB(t w ,t) only depends 
on the the time difference r = t — t w , due to the time translational invariance. On the other 
hand, in systems presenting aging, it depends on both r and t w . It is important to stress 
that the aging phenomenon is more than just the loss of the time translational invariance; it 
consists in the relaxation of Cab(£™, t w + T ) slowing down as the "waiting time" t w increases. 

The simplest aging phenomenon occurs when the two-time correlation function depends 
only on the time ratio r/t w , and it is sometimes called "full aging". This behavior is 
exhibited by some simple models, such as mean field model of spin glasses 

as 

and the 

one-dimensional Ising model at zero temperature 0, 0]. More complicated dependencies 
like lni/lntu, have also been found |6| and, with more generality, behaviors of the form 
h(t)/h(t w ), with different functions h p], Q]. 

Granular media are inherently non-equilibrium systems, due to the dissipative character 
of the interactions between grains. There is a continuous loss of kinetic energy and the system 
tends to a rest, unless energy is being continuously injected into the system, for instance 
through a vibrating wall. A kind of typical experiments carried out in dense granular systems 
are those designed to investigate compaction |8( . Usually, the system is submitted to a series 
of separated pulses or taps of a given short duration. After each tap, the system is allowed to 
relax freely until it reaches a metastable configuration with all the particles at rest. In this 
state, there are many permanent contacts between particles. Next, the system is tapped 
again and the process is repeated many times. Properties of interest like the volume or 
the energy are measured at each rest configuration. In this way, the evolution of these 
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properties as a function of time, measured in number of taps, is obtained. For large number 
of taps, the system tends to a steady state with a density that is a monotonically decreasing 
function of the tapping intensity [8|, [9|. The relaxation of the system towards the steady 
state is very slow and clearly non-exponential. Furthermore, when the intensity of tapping 
is changed cyclically, hysteresis phenomena in the density are observed . These behaviors 
are similar to those found in structural glasses when submitted to cyclical variations of their 



temperature 



10| and are a quite strong evidence of the presence of aging phenomena in 



compact granular media. Actually, aging has been observed both experimentally 11] and 
also identified in simple models of compaction I2I, Q- 

A completely different regime of granular systems are the so-called granular gases, in 
which there are not permanent contacts between particles, but they move freely and indepen- 
dently between collisions. The simplest possible state of a granular gas is the homogeneous 
cooling state (HCS), whose temperature decays monotonically in time. At a microscopic 
level, the HCS is assumed to be characterized by a phase space probability distribution in 
which all the time dependence occurs through the temperature. For the case of hard par- 
ticles, this implies a scaling property and the possibility of identifying some features of the 
time dependence of many relevant properties of the system, without carrying out explicit 
calculations. Among these properties are ensemble averages as well as time-correlation func- 



tions 



14 



151 ] . The above peculiarities render the HCS a good candidate to investigate in 



detail and by means of analytical methods the possible existence of aging and, in the case 
of a positive answer, its origin and properties. In spite of the above, it has not been until 
very recently that attention has been devoted to this particular aspect of the HCS 16] . The 
aim of this paper is twofold. Firstly, to investigate in detail the possible existence of the 
aging phenomenon in the HCS of a granular fluid of inelastic hard particles, deepening into 
its origin. This is done both analytically and by means of particle simulation methods. The 
theoretical analysis will be based on the existence of the HCS at the level of the iV-particle 
distribution function and its scaling property, for arbitrary density and inelasticity. Then, 
the accuracy of the predictions following from the theory provides strong support for the 
existence of the HCS in the context of a many body theory. This is precisely the second aim 
of the paper. 

The presentation here proceeds as follows. In Sec. [Til some general properties of the 
iV-particle distribution function defining the assumed HCS of a system of inelastic hard 
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spheres or disks are shortly reviewed. By means of an appropriate scaling of the dynamics 
of the system, general features of the time dependence of the time-correlation functions can 
be identified. This is discussed in Sec. 1111} and the results presented hold, in principle, for a 
general pair of dynamical variables and arbitrary density and inelasticity. It is shown that 
the system exhibits full aging as a direct consequence of the scaling property of the HCS. 
The particular case of the total energy of the system is considered in Sec. IIVI The time self- 
correlation function of this property in the HCS is known in detail quite accurately in the 
low density limit, then allowing to get very detailed information about its time behavior. 
Comparison of the theoretical predictions with molecular dynamics simulation results is 
also presented. A quite good agreement is observed. Finally, Sec. [V] contains some general 
comments and final remarks. 



II. THE HOMOGENEOUS COOLING STATE OF A GRANULAR FLUID 

Consider a system on N inelastic hard spheres (d = 3) or disks (d = 2) of mass m and 
diameter a. The position and velocity of particle r at time t will be denoted by q r (t) and 
v r (t), respectively. The dynamics of the system consists of free streaming, i.e. straight line 
motion along the direction of the velocity until a pair of particles, r and s, is at contact, at 
which time their velocities v r , v s change instantaneously to v' r , v' s according to 

V r = V r Y~ ' Vrs > °"' t 1 ) 

1 + a 

v' s = v s + — y~ ' Vrs ) °> ( 2 ) 
where v rs = v r — v s is the relative velocity and a is a unit vector along q rs = q r — q s at 
contact. Finally, a is the coefficient of normal restitution, defined in the range < a < 1 
and that will be considered here as a velocity- independent constant. The sequence of free 
streaming and binary collisions determines a unique trajectory of the system. 

A macroscopic state is specified, at the statistical mechanics level, in terms of a proba- 
bility density p(T,t), with T denoting a point in the 2Nd dimensional phase space of the 
system, T = {q 1: v±, . . . , q^, v N }. The macroscopic variables of interest are the average of 
microscopic observables A(T) at a given time t, defined in the two equivalent forms 

(A(t)) = j dTp(T)e tL A(T) = j dT A(T)e~ tZ ~p(T). (3) 



In the above expressions, L is the generator of the dynamics for phase functions, while L is 
the generator of the dynamics for distribution functions. Their expressions are 

N „ N N 

L ( r ) = 5>-^- +2EE T (V). (4) 

r=l r r=l s^r 

N „ N N 

r=l r r=l s^r 

with the binary collisions operators T(r, s) and T(r, s) defined by 

T{r,s) = a d - 1 J dae{-a-v rs )\a-v rs \6{q rs -(r){b rs -l), (6) 

T(r,s) = a d ~ x J dae(a ■ v rs )\* ■ v rs \ [a" 2 5(qr rs - <r)b;} - 5(q rs + tr)] . (7) 

In these expressions, da is the solid angle element corresponding to a, cr = ctct, qr rs = q r — q s , 
and B(x) is the Heaviside step function. Moreover, 6 rs is the substitution operator that 
replaces the velocities v r and v s to its right by their "postcollisional" values accordingly 
with Eqs. (CQ) and (J2]). Thus for an arbitrary function F, 

b rs F(v r ,v s ) = F{v' r ,v' s ). (8) 

Finally, the operator b~} is the inverse of b rs , i.e. it changes the velocities v r , v s by their 
"precollisional" values, 

b- 1 F(v r ,v s )=F(v';,v':i (9) 
1 + a 

V r = V r 7; " V rs) (10) 

lot 

v" = v r + (a ■ v rs ) a. (11) 

In summary, the dynamics of the probability distribution function in phase space is 
governed by the Liouville equation 

i+LW,t) = 0. (12) 

Due to the energy dissipation in collisions, there is no stationary solution to the above Liou- 
ville equation, except in the elastic limit a = 1. A granular temperature T is usually defined 
from the average of the energy density. For a homogeneous state of a system composed of 
hard particles, it is given by 

T(t) = J^W)), (13) 



with E being the total (kinetic) energy of the system. By using Eq. Q, it is found 

dT{t) 



dt 

where the "cooling rate" ((t) is identified as 



<W(t). (M) 



^=-T(m {Lm) - - (15) 

Of course, there is a large class of time-dependent homogeneous states, depending on the 
initial preparation. Here, it will be assumed that, after a few collisions per particle, there is 
a relaxation of the velocity distribution towards a "universal" form, characterized because 
its entire time dependence occurs through the cooling temperature. This special state is 
called the homogeneous cooling state (HCS) and, at the macroscopic level, it is defined 
by a uniform number density rih, a uniform but time-dependent temperature Th(t), and a 
vanishing flow velocity. Because of the absence of any additional microscopic energy scale 
for hard particles, its distribution function has the form 

Ph (T, t) = [£v (t))- Nd p* h ^, r, a = 1, . . . , tf}) , (16) 

where Vo(t) = (2T/ l /m) 1 ^ 2 is a thermal velocity and £ = (n^c ^ 1 ) -1 a characteristic length 
proportional to the mean free path. The above special form of the iV-particle distribu- 
tion function allows to determine the temperature (and time) dependence of many average 
properties without explicit calculations. This fact will be actually exploited in the following. 
The existence of the HCS solution to the Liouville equation has already been assumed 



several times in the literature 



13, 



18| . Although there is no direct proof of it, nor a con- 
structive solution of the Liouville equation for this state has been developed, molecular 
dynamics (MD) simulations have shown that some of its implications, e.g. the scaling law 
for the temperature mentioned below, are observed in detail. Additional support has been 
provided by means of a time scale change that transforms the assumed HCS distribution 



into a time-independent distribution 



201 ]. MD simulations seem to confirm the existence 
of the steady state, that is reached after a few collisions per particle. A more demanding 
evidence of the existence of the HCS with a distribution function having the scaling property 
given in Eq. (TIB]) is provided by the results to be reported in this paper. 

The temperature dependence of the cooling rate for the HCS, (h(t), can be determined 
by dimensional analysis to be (h [n h ,T h (t)) oc T h (t) 1 / 2 . Now, Eq. flU]), particularized for the 



HCS of a system of hard spheres or disks, can be integrated, to obtain the time dependence 
of the temperature 



T h {t) = T h (t') 



l Cvo(t')(t-t') 



2(: 



with 



v (t) 



(17) 



(18) 



being a dimensionless time-independent cooling rate. This algebraic decay of the tempera- 
ture of the HCS is known as the Haff law \2l\ . 

For the analysis of the HCS, it is useful to introduce the dimensionless time scale s defined 
through 



s(t) 



(19) 



o 



Therefore, s is proportional to the accumulated average number of collisions per particle in 
the time interval (0,t). In terms of this new time variable, the cooling law (FlTl) becomes 

T(s) = T(s')e-( S - S 'K* . 



The t and s time scales are related through 

2 , 
s = — In 



(20) 



(21) 



as can be directly seen by comparison of Eqs. ( fTTI) and ( 1201) or, equivalently, by direct 
integration of Eq. ffT9l . 



III. TIME-CORRELATION FUNCTIONS IN THE HCS 

As indicated in the previous section, the scaling property of the distribution function of 
the HCS implies that the time-dependence of many macroscopic properties of the system 
can be identified without carrying out explicit calculations. Let A(T) be a homogeneous 
function of degree a of the velocities of the particles. Then, it is 

A(T) = A ({q r , v r ;r = 1, . . . , N}) = A ({£q* r , v (t)v* r ; r = 1, . . . , N}) 

= v%(t)A({£q* r ,v* r ;r = l,...,N}), (22) 

where q* = q r /£ and v* = v r /v (t). Examples of this kind of properties are the center of 
mass velocity or the total energy of the system. The average value of A in the HCS is 

(A(t)) h = J dTA(T) Ph (T,t) = v a (t)(A)l (23) 
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where 

(A)* h = J dT*p* h (T*)A ({£q* r , <; r = 1, . . . , N}) , (24) 
and T* = {q*, v*;r — 1, . . . , N}. Thus all the time dependence of (A(t))h is in the factor 

Suppose next that B(T) is also a homogeneous function of the velocities of degree b, 

B (r) = v b (t)B ({£q* r , <; r = 1, . . . , iV}) , (25) 
and consider the HCS time-correlation function for A and B defined as 



C AB (t,t') = (A(t)B(t')) h - (A(t)) h (B(t')) h , 



(26) 



for t > t' > 0. By carrying out the transformation to dimensionless variables, it can be 
shown that [14] 

(A(t)B(t')) h = v a (t)v b (t'){A(s - s')B)t (27) 

with 

(A(s)B)* h = J dT* p* h (T*)A ({e q ;„ <} , a) B ({£q* r , <}) . (28) 

Here 

A ({lq* r , <} , 8 ) = e sC *A ({£q* r , <}) , (29) 
where £* is the new generator for the dynamics of the phase functions, 



N 



r=l 



d 

dv* 



L*(T*), 



L*(T*) 



v (t) 



^( r ) = [L(r)] {qr=q * >Vr=v * }>a=a * 



(30) 



(31) 



The first term on the right hand side of Eq. (130]) . is due to the time- dependent scaling of 
the velocities with v (t). Use of Eqs. (I2"3"j) and (f2"T|) into Eq. fTSBI) yields 

C AB (t,t') = v a (t)v b (t')C* AB (s - 0, (32) 

C AB (s) = (A(s)B)l-(Ay h (B)l. (33) 

To identify the aging phenomena clearer, it is convenient to normalize the correlation func- 
tion to unity for t = t', by defining a relaxation function (f)AB{t,t') as 



4>AB(t,t') 



C AB (t',t>) 



vo(t) 
vo(t') 



'AB 



s-s'-C 



- {s - s,) ^<Pab(s-s';C). 



(34) 



Upon writing the last equality above, use has been made of Eq. (1201) . The dimensionless 



relaxation function <p* AB above is 



_ (A(s)B}l-(A}l(B} h 



It follows from Eq. (jMl) that the two-time correlation function depends on time only through 
the difference s — s'. Consequently, there is no aging when time is measured in the dimen- 



sionless scale s. In ref. [la], the aging property of the velocity time-autocorrelation function 
of a granular gas of inelastic hard particles was investigated by means of MD simulations. 
Time was measured by the average cumulated number of collisions per particles that, as said 
above, is proportional to the dimensionless time scale s used here. The simulations indicate 
that the velocity autocorrelation function, C vv (s, s') in the language used here, depends both 
on s' and s—s', and from this feature the authors conclude that the system exhibits aging. Of 
course, the analysis developed here applies for the case of the velocity time-autocorrelation 
function, corresponding to a = b = 1 and, therefore, aging should not be expected according 
to the results derived above. This apparent discrepancy seems to occur because of a wrong 
use of the physical concept of aging in ref. [if]]. As pointed out in the Introduction, for the 
existence of aging, it is not enough the dependence of the time-correlation function on both 
s' and s — s'. In fact, the analysis developed here leads to 

Cab(s, s') = C AB (s', s')4> AB (s, s') = f{s')e-^ s - s '^\ B {s - s'; (*). (36) 

Nevertheless, all the dependence on s' occurs in the prefactor f(s') = Cab{s' ', s') and, there- 
fore, there is no real aging, since the decaying rate is always the same and only the initial 
value changes with s' for constant s — s'. 

On the other hand, the aging phenomenon shows up in the original time scale t in the 
limits Ch(0)t' ^> 1 and (h(0)t ^> 1. In this regime, it follows from Eq. ( f2Tl) that 

s-s'~llni, (37) 
and, therefore, Eq. (13^1) takes the form 



4>AB(t,t') 



'AB 



2 , t 

^ln-;C 
C* t 



= F[^C), (38) 



y 

valid for t, t' 3> Gi(0) -1 - This result is the mathematical expression of the aging phenomenon. 
The normalized time-correlation function depends on the initial and final times, t' and t, only 



through their quotient t/tf, so the system exhibits full aging, as defined in the Introduction. 
The remaining parameter determining the long time behavior of the correlation function 
is the dimensionless cooling rate (*. It is worth to stress the generality of this result. 
No limitation on the degree of inelasticity or density has been introduced. In fact, both 
magnitudes are relevant in determining, through the value of C(0), the time region in which 
the aging behavior predicted by Eq. (I36I) is to be expected. The only hypothesis made 
here is the existence of the HCS with a distribution function having the scaling form given 
in Eq. (fl6|) . Of course, the system is assumed to stay in that state for all the relaxation 
time t considered. This requires that the HCS, in addition to exist, be stable, at least in 
a determined region of parameters. The extensive measurements of the velocity correlation 



221 ] show that it is possible to cover 



function by means of MD simulations reported in ref. 
a wide range of density and inelasticity in which the observed homogenous state appears 
to be stable. This is confirmed by the simulation results to be presented here in the next 
section. 



IV. ENERGY TIME CORRELATION FUNCTION IN THE HCS FOR A DILUTE 
GRANULAR GAS 

Theoretical predictions for the explicit form of the function 0^g(s,C*) are scarce in the 
literature. An exception is the total energy autocorrelation function, C EE (t, t'), for a dilute 
granular gas of hard particles. By projecting the Liouville equation onto the hydrodynamic 



modes in the low density limit, an analytical expression for C EE was derived in ref. [15| 



C EE (t, = NT(t)T(t')e(a)e^ s - s '^ , (39) 

where e(a) is a given function of only the coefficient of restitution. Then, 

, /. ,,s _ C EE (t,f) T(t) _( fl - B ')£ / 4n x 

Since the energy E is a homogenous function of degree a = 2 of the velocities of the particles, 

comparison of the above expression with the general result given in Eq. (|34p leads to the 
identification 

4> EE (s;C) = e-^ s . (41) 
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Therefore, using Eq. (I38I) . it follows that the aging behavior of the total energy of a dilute 
granular gas of hard particles is given by 

<WM')=e-^ (s - s,) ~ (Q 3 , (42) 
for t,t' ^> (^(O) -1 . An explicit expression for the cooling rate of a dilute granular gas has 



been derived from the Boltzmann equation in the so-called first Sonine approximation 23]. 



results obtained by 



241 ] in the thermal 



This expression has been shown to accurately agree with the numerica 
means of the direct simulation method of the Boltzmann equation 2C 
velocity region, i.e. for velocities of the order of vo(t). 

Due to the simplicity of <fr EE in the present case, the asymptotic form of the normalized 
time correlation function for the energy, <pEE(t,t'), does not depend on the value of £*, 
what implies that it is independent of the density rih and of the coefficient of restitution a. 
Consequently, if <pEE(t, t') is plotted as a function oit/t' for times in the range t, t' 3> C(0)~ 1 , 
the curves corresponding to different densities and inelasticities should tend to collapse on 
a unique one. Of course, the density range to which this result applies is restricted by the 
validity of the Boltzmann description. 

To check the above theoretical predictions, we have performed MD simulations of a system 
of inelastic hard disks (d = 2). In Fig. [I], results obtained for a system of N = 10 3 particles 
with coefficient of normal restitution a = 0.95 and density rihcr 2 = 0.02 are reported. For 
these values of the parameters, the HCS is stable, since the critical length is larger than the 
size of the system and, therefore, the velocity vortices and high density clusters characteristic 
of the clustering instability [25] cannot develop. Of course, in all the simulations it has been 
verified that the system remains homogeneous. The reported results have been averaged 
over 1200 trajectories of the system. The several curves correspond to different values of 
0i(0)t', as indicated in the figure. The value of C/i(0) has been estimated by using the low 



density expressions derived in [23 1. 



In agreement with the analysis presented here, it is observed that as t' increases the time- 
correlation function approaches a form that depends only on the value of the time ratio 
t/t'. More precisely, for 0,(0)*' > 24, all the plotted curves coincide within the statistical 
uncertainties. Note that the latter increase as the value of t grows, for a given value of t'. 
Moreover, the asymptotic curve agrees with the theoretical prediction given by Eq. (|42l) . 
whose graphical representation is the solid line in the figure. This good accuracy is consistent 

11 




(t'/t) 



FIG. 1: Dimensionless time self-correlation function of the total energy, <f>EE(t,t'), for a dilute 
granular gas of inelastic hard disks in the HCS. The density is n^c 2 = 0.02 and the restitution 
coefficient a = 0.95. The different symbols correspond to different values of the initial time, as 
indicated. The solid line is the theoretical prediction describing the full aging phenomenon, Eq. 
121). 



with the value of the density of the system, that is small enough as to expect a low density 
description, at the level of the Boltzmann equation, to apply. 

In Fig. [2J a similar plot is given, but now two systems, one with n^cr 2 = 0.02, a = 0.85, 
and the other with rihCT 2 = 0.1, a = 0.98, are considered. For the sake of clearness, only 
simulation data corresponding to large waiting times have been included. Again, a good 
agreement with the behavior predicted by Eq. (H21 is observed. The same behavior has been 
obtained for other densities between the two above values. This confirms the independence 
of 4>EE{t, t') from the density and the inelasticity at low density. In fact, the good agreement 
observed for rihcr 2 = 0.1 must be stressed, because at this density the gas can not be 
considered as very dilute. 

Prompted by the results for r^a 2 = 0.1 given in Fig. [2], the correlation function <f>EE(t,t') 
has also been evaluated at a definitely non small density, namely n^u 2 = 0.2. For this 
value, density corrections to the low density behavior are clearly identified in most of the 
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FIG. 2: The same as in Fig. [H but for different values of the (low) density and of the restitution 
coefficient. Filled symbols correspond to a system with n^cr 2 = 0.02 and a = 0.85, while empty 
symbols refer to a system with n^cr 2 = 0.1 and a = 0.98. 



equilibrium properties of a molecular gas. It must be mentioned that, in order to keep 
the system well inside the stable region of parameters and with a number large enough of 
particles, the value of a must be rather close to unity. Moreover, they have been averaged 
over 1500 trajectories. The results shown in Fig. [3] have been obtained with a = 0.98 
and N = 700. Once again, a tendency towards a behavior depending only on the ratio 
t/t' as t' increases is clearly identified, i.e. the system exhibits full aging. Besides, and 
rather surprisingly, the aging phenomenon seems to be accurately described over several 
decades by the law (t'/t) 3 , that was obtained here in the context of very dilute granular 
gases (Boltzmann limit). 



V. DISCUSSION 



The objective here has been to explore the existence of aging in the homogeneous cooling 
state of a granular gas. This has been done by exploiting the assumed scaling property of 
the iV-particle distribution of this state for a system of inelastic hard spheres or disks. In 
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(t'/t) 



FIG. 3: The same as in Fig. [TJ but for n^u 2 = 0.2 and a = 0.98. A clear tendency towards a 
function depending only on t/t', characteristic of full aging, is observed as t' increases. 

fact, the presence of aging and its specific form turn out to be directly associated to the 
scaling property of the distribution Q] • From this perspective, the agreement between 
molecular dynamics simulations and the theory discussed in this paper, is an almost direct 
proof of the existence of the homogeneous cooling state also at the level of the full many 
body pseudo-Liouville equation. This would extend the rather well established fact that 
the inelastic Boltzmann and Enskog equations have such a solution, consistently with some 



previous results 



19 



20J. 



Granular media are inherently non-equilibrium systems due to the lack of energy conser- 
vation in the interactions between grains. They present a very rich phenomenology which, 
sometimes, is similar to that of normal, molecular systems. Moreover, it has been verified in 
the last years that the methods of kinetic theory and non-equilibrium statistical mechanics 
developed for normal fluids, can be extended to granular fluids, yielding to results having an 
analogous structure. There are also significant differences, but they are well understood as 
consequences of the inelasticity. One characteristic feature of granular systems is that quite 
often they exhibit the phenomena in a much simpler context than molecular systems. This 
refers to both, theoretical and experimental views. In this paper, the simplicity of the HCS 
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of a granular fluid has allowed the identification of aging and the derivation of its explicit 
form in the dilute limit, Eq. (|42p . for the correlation of the total energy of the system. The 
analytical expression, that corresponds to the so-called full aging, has been shown to be 
in perfect agreement with MD simulation results. In particular, the long time limit of the 
normalized time-correlation function of the total energy is independent of the inelasticity 
and the density. Quite interestingly, the simulations indicate that this independence seems 
to extend to densities beyond the dilute limit. 

To really appreciate the results presented here, it is important to differentiate between 
both the existence of aging and the specific, particular, form of the law governing it. Equa- 
tion f!38p implies the existence of full aging in the system, i.e. the normalized time-correlation 
function of the properties A and B is not a function of the time difference t — if and depends 
on the time ratio t/t' . The only necessary condition to derive this equation is the existence 
of the HCS itself, as discussed above. On the other hand, identification of the function F 
in Eq. ( |38l) requires more detailed additional analysis, which has been carried out only in 
the low density limit up to now. The above leads to Eq. ( H2l) in the particular case of 
the properties A and B being both the total (kinetic) energy of the system. In this sense, 
the results presented in Fig. [3] strongly support the existence of the HCS and the scaling 
property of its iV-particle distribution function at high densities, independently of whether 
or not the convergence occurs towards the power law given by Eq. fj42|) . as suggested by the 
simulation results. 

In real granular gases, the restitution coefficient is not constant, but it depends on the 
impact relative velocity. Then, the distribution function of the HCS does not scale in the 
form given by Eq. ([16]) and, therefore, the discussion in the present paper does not apply in 
principle, although it can provide an accurate approximation to the actual behavior of the 
system. To be more precise, consider a given model of granular gas with velocity-dependent 
coefficient of normal restitution. Now, the problem being addressed has two energy scales. 
One is the total energy per particle or, equivalently, the cooling temperature Th(t). The 
other energy scale, e, is fixed by some property of the specific collision model. Define a 
dimensionless parameter 

For hard spheres, e = and so e* = 0. It is in this limit when the distribution function of 
the HCS has the scaling property (jTEj) [261 ] . For e* > 0, the scaling does not hold exactly, 
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but it can be an appropriate description as long as e* <C 1, i.e. the interaction be sufficiently 
hard and/or the kinetic energy of the particles be sufficiently large. 
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